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Methane, the primary constituent of natural gas, binds too weakly to nanostruc- 
tured carbons to meet the targets set for on-board vehicular storage to be viable. 
We show, using density functional theory calculations, that replacing graphene by 
graphene oxide increases the adsorption energy of methane by 50%. This enhance¬ 
ment is sufficient to achieve the optimal binding strength. In order to gain insight 
into the sources of this increased binding, that could also be used to formulate de¬ 
sign principles for novel storage materials, we consider a sequence of model systems, 
that progressively take us from graphene to graphene oxide. A careful analysis of 
the various contributions to the weak binding between the methane molecule and 
the graphene oxide shows that the enhancement has important contributions from 
London dispersion interactions as well as Debye interactions and higher-order electro¬ 
static interactions, aided by geometric curvature induced primarily by the presence 
of epoxy groups. 
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I. INTRODUCTION 


Though most vehicles today run on gasoline, there are several pressing reasons to switch 
instead at least partly to natural gas - it is cheaper as well as cleaner burning, and reserves 
of natural gas are distributed more uniformly across the globe. At present, vehicles running 
on natural gas make use of CNG (compressed natural gas) technology, involving the use of 
heavy, unwieldy and dangerous gas cylinders, and high pressures. The hope is to replace 
this by adsorptive storage in a solid state material. However, achieving the target storage 
capacities and binding strengths that allow one to compete with gasoline-based technologies 
poses a formidable materials and engineering challenge.^ 

High storage capacities for methane (the primary constituent of natural gas) have been 
achieved with metal organic frameworks P® however, at present these materials are rather 
expensive to be used at an industrial scale, and it is worth exploring other classes of materi¬ 
als in parallel. While nanostructured carbons (such as graphene and carbon nanotubes) are 
also attractive candidates for use as adsorptive natural gas (ANG) sorbents,^ the binding 
of methane to these systems is too weak for on-board vehicular applications The exper¬ 
imentally measured value of E^ds, the adsorption energy of methane, is ~12 kJ/mol on 
graphite,!^ which is well below the estimated optimal value of 18.8 kJ/mol.^ Calculated val¬ 
ues for A'ads on graphene range from 11-16 kJ/mol, depending on the method usedNote 
that the target value for E^^ds falls between the binding strengths typical of physisorption 
and chemisorption. 

In addition to E^ds, other parameters that determine gas storage capacities are surface 
area, void space and density.l^ Previous authors have explored the possibility of increasing 
the uptake of methane in carbon-based systems by, e.g., tuning pore size,^ and increasing 
porosity by pyrolysis.l^^ Here, we will focus however on just one parameter, viz., A'ads- Di¬ 
verse strategies have been attempted to tweak A'ads upward - by e.g., edge-functionalizing 
by various chemical groups,! ^ * ^^ * ^^ * activation by steam and carbon dioxide,^ and the in¬ 
troduction of defects and curvature.!^ However, hitherto, none of these have succeeded in 
enhancing A'ads sufficiently. Here, we show, using ab initio density functional theory calcula¬ 
tions, that replacing graphene by graphene oxide succeeds in increasing A'ads sufficiently so 
as to meet the establishecP target value. We also show that this enhancement has important 
contributions from three different effects, alf of which act in concert here. We note also that 


2 




graphene oxide has the additional advantage that it is more easily synthesized by a variety 
of chemical routes, compared to graphene .1221 

There has been earlier theoretical work on hydrogen binding on graphene oxide (GO)j22E31 
however, these studies used the GO as a substrate for metal adatoms such as Ti and Ni, which 
act as the binding centers for gas molecules; we note that adding such metal atoms raises 
the weight of the system. Experimentally synthesized graphene oxide based framework^ 
show high isoteric heat of adsorption and enhanced uptake of H 2 , due to the large surface 
areas of frameworks achieved by boronic acid based pillaring between GO layers. 


II. CALCULATION DETAILS 

Our ab initio density functional theory (DFT) calculations have been performed using 
the Quantum ESPRESSO package,l2^ with plane wave cut-offs of 40 Ry and 400 Ry for the 
wavefunction and charge density, respectively, and ultrasoft pseudopotentials .^221 Exchange 
correlation interactions were treated within a generalized gradient approximation .^221 Im¬ 
portantly, in order to have an accurate treatment of weak binding at a relatively modest 
computational cost, we have incorporated van der Waals interactions using the “DET-D2” 
method.®' The Brillouin zone was sampled with grids commensurate with a (12 x 12 x 1) 
mesh for the primitive unit cell of graphene. Convergence was aided by using cold smearin^^Ql 
with a width of 0.001 Ry. Artificially periodic images normal to the plane of the graphene 
sheet were separated by 20 A. All atomic coordinates were allowed to relax, using Hellmann- 
Feynman forces, with a convergence threshold of 10“^ Ry/bohr. Several initial orientations 
of the methane molecule relative to the substrate were considered, in order to span the space 
of possible geometries. 

The adsorption energy of methane is given by: 

A'ads (A'sub-|-CH4 A'sub ^^0114); (1) 

where E'sub+CH 4 , A'sub and E'ch 4 are the calculated total energies of the substrate with ad¬ 
sorbed methane, the substrate alone, and an isolated methane molecule in the gas phase, 
respectively. 
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III. RESULTS 


First, we consider methane adsorption on bare graphene; we label this confignration C(I). 
We obtain E^ds = 14.37 kJ/mol; this falls within the range of previons theoretical valnes.^^^ 
Next, we wish to consider the adsorption of methane on graphene oxide. The structure 
of graphene oxide continues to be a matter of debate, and several models exist for it. The 
carbon atoms are arranged in the honeycomb lattice characteristic of graphene. However, 
hydroxyl (-OH) and epoxy (-0-) groups are not attached in a periodic manner to these 
carbon atoms, so that the overall structure lacks crystalline symmetryj ^^F2 | 

We first wish to examine the energetics of methane adsorption on areas of GO that are 
densely covered with functional groups. To do this, we make use of the periodic graphene 
oxide (PGO) approximant of Wang et al, which features a repeating motif of a triplet of 
functional groups (two -OH on one side of the graphene sheet, and one -O- on the other 
side) attached to a sixfold carbon ring; this model was motivated by the results of nuclear 
magnetic resonance experiments,^^ and guided by the energetics of various configurations 
as obtained from DFT. The optimal geometry obtained when a single methane molecule 
is adsorbed within the (2\/3 x \/3) primitive unit cell of PGO is shown in Figs. (a) and 
(b); it leads to Eads = 17.04 kJ/mol. In the more dilute case, where a single methane 
molecule is adsorbed instead in a (2^/3 x 2\/3) cell [see Figs. (c) and (d)], Eads decreases 
to 15.71 kJ/mol, due to a lessening of the attractive interaction between methane molecules 
in neighboring cells. 

While these values of Eads do not quite meet the target^ considered desirable for on-board 
storage, we must remember that the PGO approximant is not a realistic model of GO, since 
it describes a substrate that is densely covered with functional groups, whereas graphene 
oxide is believed, in reality, to contain in addition, considerable amounts of patches of bare 
graphene, with many interfaces present between covered and bare areasAccordingly, we 
next consider adsorption of methane at an interface between PGO and bare graphene. We 
model this interface as either a zigzag or armchair-like line. -OH groups near the interface 
now have considerable orientational freedom, and are in fact found to re-orient themselves, 
compared to their alignments in the PGO structure. 

For a methane molecule adsorbed in a (4\/3 x 6) cell containing a zigzag interface, we 
find two stable adsorption geometries that are very close in energy, we refer to these as 
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(a) ( 2 V 3 X V3) PGO (b) ( 2 V 3 X ^3) (c) (2^3 x 2^3) PGO (d) (2^3 x 2^3) 

top-view PGO side-view top-view PGO side-view 



(e) Z-INNER top-view 



(f) Z-INNER side-view 


(g) Z-OUTER 


top-view 



(h) Z-OUTER side-view (i) A-OUTER Top 


(j) A-OUTER side-view 


view 


FIG. 1. Optimal relaxed geometries and charge redistribution Ap for methane adsorbed on 
(a)-(d) PGO and (e)-(j) at an interface between bare graphene and PGO. Isosurfaces of Ap = 
±2.8 x 10 ^ e/bohr are plotted, Red and green lobes denote the gain and depletion of electronic 
charge, respectively. Color scheme for atomic spheres: H (black), C (gray), O (blue). 


“Z-INNER” [see Figs. ^ (e) and (f)] and “Z-OUTER” [see Figs. [T| (g) and (h)]; in these 
two configurations the methane molecule is positioned near the interface, but over covered 
and bare regions of graphene, respectively. Most interestingly, adsorption near the interface 
increases the binding, pushing Eads up to 21.70 kJ/mol and 21.59 kJ/mol for the two geome¬ 
tries, respectively. Similarly, for adsorption near an armchair interface in a {2^/3 x 4\/3) cell 
[see Figs.[^(i) and (j), labeled “A-OUTER”], the value of Eads is again increased, this time 
to 20.87 kJ/mol. These three values of Eads uU meet the target set for on-board storage to 
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(a) C(I) side-view 


(b) C(II) side-view (c) C(III) side-view 


(d) C(IIIa) 
side-view 



side-view 


FIG. 2. Optimal relaxed geometries and charge redistribution Ap for methane adsorbed on 
bare graphene [configuration C(I)] and a sequence of model systems [configurations C(II)-C(VI)]. 
Isosurfaces of Ap = ±2.8 x 10“^ e/bohr are plotted, Red and green lobes denote the gain and 
depletion of electronic charge, respectively. Color scheme for atomic spheres: H (black), C (gray), 
O (blue). 


be viable,® suggesting that graphene oxide may be a good candidate material for adsorptive 
storage of natural gas. Note also that by going from graphene to graphene oxide, we have 
raised the binding strength by one and a half times. 

We now proceed to examine the origins of this enhanced binding, and try to break it down 
into various energetic and structural components. The main energetic contributions relevant 
to weak binding such as that being considered here, are London dispersion interactions 
(LDI), and various electrostatic interactions (some of which are subsumed into van der 
Waals interactions, such as the Debye interaction between permanent and induced electric 
dipoles). The part of E^, arising from LDI is given by By = - -BoS - 

where the three terms on the right-hand-side are the contributions from LDI to the total 
energies of the corresponding systems. We find that the values of for Z-INNER, Z- 
OUTER and A-OUTER are 19.08, 21.58 and 19.05 kJ/mol, respectively. This shows that 
the majority of the binding comes from LDI, and thus any theoretical treatment that does 
not account properly for these may be expected to be inaccurate. However, there are also 
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some repulsive contributions to -Eads, and the LDI is not the only attractive interaction that 
offsets these. In fact, we will see that attractive electrostatic interactions between induced 
charges on the methane molecule and the dipolar -OH and -O- groups provide important 
contributions toward stabilizing the adsorbed methane, and toward pushing Eads into the 
desired window. As an illustration that interactions other than LDI are important, we note 
that on going from CH 4 adsorbed on PGO in the {2^/3 x 2 \/ 3 ) cell, to the Z-INNER interface 
configuration, Eads goes up by 5.99 kJ/mol, whereas El is increased only by 1.37 kJ/mol. 

In Fig. in addition to depicting the relaxed structures of the methane-|-substrate sys¬ 
tems, we have also plotted isosurfaces of the charge redistribution Ap = Psub+CH 4 — PCH 4 — 
Psub, where p is the electronic charge density, and all three terms on the right-hand-side are 
evaluated at the relaxed geometry of the combined system. Red and green lobes correspond 
to an accumulation and depletion of electronic charge, respectively. In Fig. we can see 
very clearly that though methane in the gas phase possesses neither a dipole nor quadrupole 
moment, when adsorbed on PGO, it displays a considerable dipole moment. This dipole 
moment is larger when it is adsorbed at an interface than on PGO; a visual indication of 
this can be obtained by observing the larger red and green lobes in Figs. (e)-(j) than in 
Figs.|3(a)-(d). 

In order to break down the interaction of the methane molecule with graphene oxide into 
various components, we now consider a sequence of simple model systems, designed so as to 
progressively take us from bare graphene to graphene oxide, while allowing us to separate 
out the various effects at play. 

We have already mentioned configuration C (I), where methane binds to a clean graphene 
sheet. In this case, the binding is almost entirely due to LDI; electrostatic interactions are 
negligible, as evidenced by the very small red and green lobes in Fig. [^a). 

Next, we consider graphene functionalized with a single epoxy group. It is interesting 
to note that it is far more favorable for CH 4 to bind on the opposite side of the epoxy 
[configuration C(III), see Fig. [^c), Eads = 17.11 kJ/mol] than on the same side [C(II), see 
Fig. ib), Eads = 12.69 kJ/mol]. This preference has a structural origin: the attachment 
of the epoxy group results in a considerable local buckling of the graphene sheet, with the 
two G atoms that are bonded to the O atom being displaced vertically by ~ 0.24 A relative 
to the rest of the atoms. It is known that methane molecules prefer to bind to a concave 
curvature of a graphene sheetpl this is the phenomenon being observed here. If we remove 
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the O atom, but freeze the graphene sheet in the buckled configuration it assumes in the 
presence of the epoxy group [C(IIIa), see Fig. [^d)], we obtain E^ds = 15.76 kJ/mol, i.e., 
we recover most of the enhanced binding. As expected, C(III) and C(IIIa) have almost 
identical values for E^. Note also that the isosurfaces in Figs. (c) and (d) do not display 
any signature of significant electrostatic contributions to binding. 

Functionalizing graphene with a single hydroxyl group leads to a quite different scenario. 
Methane prefers to bind on the same side as the -OH [configuration C(IV), Figs, [^e)], and 
there are noticeable induced moments on the molecule, as well as some redistribution of 
charge on the substrate. The value of E^ds goes up now to 20.56 kJ/mol. Interestingly the 
induced dipole on methane is oriented oppositely, with respect to the functional group, for 
epoxy and hydroxyl [compare Figs, [^b) and[^e)]. 

Finally, we consider what happens when both an epoxy and a hydroxyl are present 
together, as in graphene oxide. In Fig. [^f) we show configuration C(V), where the -O- and 
-OH are adsorbed on opposite sides of the graphene sheet in a very low energy geometry.l^ 
The value obtained for F'ads is now increased further, to 21.78 kJ/mol. In the hypothetical 
configuration C(Va), we retain the -OH group of C(V), but remove the -0-; however, we 
freeze the system geometry at that of C(V) [see Fig. [^g)]. F'ads is almost unchanged, being 
equal to 21.66 kJ/mol, again confirming that the main role played by the epoxy group in 
enhancing binding is through the route of inducing a buckling of the graphene sheet. We 
also consider C(VI) [see Fig. [^h)], which contains both -O- and -OH, but on the same side 
of the sheet. While this is a higher energy configuration than that depicted in C(V), it leads 
to the highest binding strength of all the configurations considered in this study, giving F^ads 
= 22.94 kJ/mol. This can be readily understood as a synergistic effect of the hydroxyl and 
epoxy groups, since the electric fields of the hydroxyl and epoxy groups add together to 
result in a larger induced electric dipole on the methane molecule. 

To summarize the results we obtain in this analysis of adsorption energies, we have plotted 
our results for E^ds and El, for the various model configurations C(I) to C(VI), in the bar 
chart in Fig. [^ Note that, unlike our results for the cases Z-OUTER, Z-INNER, and A- 
OUTER, in all the cases considered here, El is larger than Eads; this is because while Eads 
contains attractive contributions from El, as well as from electrostatic interactions (which we 
will call Ees), there are also repulsive contributions. A small height difference between Eads 
and El means that Ees is large, as can be readily verified by comparing the heights of the 


bars in Fig. [^with the sizes of the red and green lobes in Figs. [T] and With the exception 
of configurations C(IIIa) and C(Va), which were hypothetical, non-realistic configurations 
introduced for the purpose of separating out the effects of buckling from the electronic or 
charge-transfer effects due to -0-, we see that as we proceed from C(II) to C(VI), Flads 
increases at every step. As we go from C(II) to C(III), Fiads increases because increases; 
this in turn is because of the buckling of the graphene sheet, caused by the epoxy group. On 
proceeding from C(III) to C(IV), El stays about the same; however, E^^g increases, because 
electrostatic attraction [which was practically absent in C(III)] increases. C(V) differs from 
C(IV) in that there is an epoxy on the bottom side (opposite the methane); this results 
in an additional buckling of the graphene sheet, as a result of which El increases, and 
thus Eg^s increases by a corresponding amount. Finally, on going from C(V) to C(VI), El 
increases slightly, but also the electrostatic interactions increase because of the synergistic 
effect already discussed above, and thus we get the largest value of Eg^. However, it is 
C(V) whose structure is both lower in energy as well as more representative of GO, and the 
value of Flads in C(V) is similar to those obtained by us for the Z-INNER, Z-OUTER and 
A-OUTER geometries. 

IV. CONCLUSION AND SUMMARY 

In conclusion, we have shown, by DET calculations, that replacing graphene by graphene 
oxide increases the strength of binding to methane by about 50%, which is sufficient to 
bring the adsorption energy into the target range for on-board vehicular storage. The 
enhancement comes from a synergy between various contributions to weak binding, viz., 
London dispersion forces and electrostatic interactions, aided by structural distortion of the 
graphene sheet by epoxy groups. This suggests that graphene oxide and graphite oxide 
should be good candidate materials for the on-board storage of natural gas. The interfaces 
between patches of densely covered GO and graphene are optimal binding sites for methane. 
Thus, if one could experimentally synthesize graphene oxide so as to deliberately engineer a 
high density of labyrinth-like interfaces between bare graphene and densely covered graphene 
oxide at the nanoscale, this should be optimal for high-capacity methane storage. Moreover, 
the understanding gained in this study can be leveraged to develop design principles that can 
be used to engineer future materials for methane storage, so as either increase the binding 
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FIG. 3. Bar chart indicating step-wise enhancements in adsorption energy F'ads) and the contribu¬ 
tion El to it from London dispersion interactions, for methane adsorbed on bare graphene (I) and 
a sequence of model systems (II)-(VI). See the text for a description of the systems. Configurations 
(Ilia) and (Va) are hypothetical in that they do not correspond to favorable configurations of the 
substrate, and are therefore shaded differently. Note that the solid blue bars increase at each stage, 
except for (II), which is energetically unfavored. 

strength in nanocarbons even further beyond the values achieved here, or to transfer this 
knowledge to other materials. 
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